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We investigate the spatial extension of weakly bound Ne and C isotopes by taking into account the 
pairing correlation with the Hartree-Fock-Bogoliubov (HFB) method and a 3-body model, respec- 
tively. We show that the odd-even staggering in the reaction cross sections of ■^'^'^I'-^^jv^g ^^^^ i4,i5,16q 
are successfully reproduced, and thus the staggering can be attributed to the pairing anti-halo effect. 
A correlation between a one-neutron separation energy and the anti-halo effect is demonstrated for 
s- and p- waves using the HFB wave functions. 
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Interaction cross sections aj and reaction cross sec- 
tions aji of unstable nuclei have been measured using 
exotic isotope beams produced by projectile fragmenta- 
tion in high-energy heavy-ion collisions. An important 
motivation for these studies is to determine the size of 
unstable nuclei [ll-Q. These measurements have in fact 
discovered the largely extended structure of unstable nu- 
clei such as "Li Q, "Be Q, and "'^^C '^]. Nuclei with 
an extended density distribution are referred to as "halo" 
nuclei and their structure has been studied extensively as 
one of the characteristic features of weakly bound nuclei. 
A proton halo has also been claimed for due to a large 
enhancement of the quadrupole moment Q. The root- 
mean-square radius diverges for s and v waves as the 
single-particle energy approaches zero [a| and the halo 
structure can be ascribed to an occupation of an / = or 
I ~ I orbit by the valence nucleons near the threshold 0] . 

Recently, Takechi et al. measured the interaction cross 
sections of neutron-rich Ne isotopes ^*~'^^Ne @ and re- 
ported a halo structure of the '^^Ne nucleus. The halo 
structure of ■^^Ne was also reported in the one-neutron 
removal reaction by Nakamura et al. Q . The interaction 
cross section aj is almost the same as the reaction cross 
section an for unstable nuclei since the cross sections for 
inelastic scattering are in general small [13, [Hi • In addi- 
tion to large reaction cross sections for neutron-rich Ne 
isotopes, the experimental data also show that the re- 
action cross sections of ^°Ne and ^^Ne are quenched as 
compared to the neighboring odd Ne isotopes, that is, a 
clear odd-even staggering has been found in the reaction 
cross sections 8']. 

In this paper, we discuss the relation between the odd- 
even staggering in the reaction cross sections and the so 
called pairing anti-halo effect [l^l- The pairing corre- 
lation has been known important in nuclei throughout 
the periodic table, giving an extra binding for paired nu- 
cleons [H [l|. It plays an essential role especially in 
loosely bound nuclei such as ^^Li and ®He (so called the 
Borromean nuclei), because the nuclei will be unbound 
without the pairing correlations. It is known that many- 
body correlations strongly modify the pure mean-field 
picture of loosely bound nuclei. It may happen that the 



strong pairing correlations do not allow the decoupling of 
weakly bound nucleons from the core nucleus and prevent 
the growth of halo structure. This pairing "anti-halo" 
effect has been suggested with the Hartree-Fock Bogoli- 
ubov (HFB) method 12[, but a clear experimental sig- 



nature of the anti-halo effect has not yet been obtained so 
far. We will examine in this paper whether the odd-even 
staggering observed in the recent experimental data for 
the reaction cross sections of neutron-rich Ne isotopes is 
an evidence for the pairing anti-halo effect. 

In order to address this question, we first study 
neutron-rich Ne isotopes using the HFB method. The 
asymptotic behavior of a Hartree-Fock wave function for 
s— wave reads 



^j{r) ~ exp(— ar). 



(1) 



where a = ^2ra\e\IW- with the HF energy e. The mean 
square radius of this wave function is then determined as 
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which will diverge in the limit of vanishing separation 
energy |e| — >■ 0. This divergence will occur not only for 
s— wave but also for p— wave, although the dependence 
on |e| is different [61 . In contrast, the upper component 
of a HFB wave function, which is relevant to the density 
distribution, behaves [l^ 



v{r) (X exp(— /3r), 



(3) 



where /? is proportional to the square root of the quasi- 
particle energy E, 
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A being the chemical potential. If we evaluate the quasi- 
particle energy in the BCS approximation or HFB with 
canonical basis, it is given as 



= ^(e - A)2 + A2, 



(5) 
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where A is the pairing gap. For a weakly bound single- 
particle state with e ~ and A ^ 0, the asymptotic be- 
havior of the wave function v(r) is therefore determined 
by the gap parameter as, 



v{r) oc exp 




(6) 



The radius of the HFB wave function will then be given 
in the limit of small separation energy |e| — as 
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As we show, the gap parameter A stays finite even in 
the zero energy limit of e with a density dependent pair- 
ing interaction. Thus the extremely large extension of a 
halo wave function in the HF field will be reduced sub- 
stantially by the pairing field and the root-mean-square 
(rms) radius of the HFB wave function will not diverge. 
This is called the anti-halo effect due to the pairing cor- 
relations [12|]. 

We now numerically carry out mean-field calculations 
with a Woods-Saxon (WS) potential and also HFB cal- 
culations using the single-particle wave functions in the 
WS potential. As examples of s-wave and p-wave states, 
we choose the 2si/2 state in ^'^O and 2p3/2 state in 
•^""^Ne, respectively. Although ^"'^Ne is most likely a de- 
formed nucleus \Wi [l7j . for simplicity we assume a spher- 
ical Woods-Saxon mean-field potential. Notice that a 
Woods-Saxon potential with a large diffuseness parame- 
ter a yields the 2p^/2 state which is lower in energy than 
the 1/7/2 state, as was shown in Ref. |18|]. We use a 
similar potential with a=0.75 fm as in Ref. [3l for '^^Ne, 
while that in Ref. 19] for ^-^O. For the HFB calculations, 
we use the density-dependent contact pairing interaction 
of surface type, in which the parameters are adjusted 
in order to reproduce the empirical neutron pairing gap 
for '^°Ne [2fl]. While we fix the Woods-Saxon potential 
for the mean- field part, the pairing potential is obtained 
self-consistently with the contact interaction. 

The top panel of Fig. 1 shows the mean square radius 
of the 2si/2 state for ^'^O, while that of Fig. 2 shows 
the mean square radius of the 2p3/2 state for ^-'^Ne. In 
order to investigate the dependence on the single-particle 
energy, we vary the depth of the Woods-Saxon wells for 
the si/2 and P3/2 states for ^'^O and "^^Ne, respectively. 
The dashed lines are obtained with the single-particle 
wave functions, while the solid lines are obtained with 
the wave function for the canonical basis in the HFB 
calculations. One can see an extremely large increase of 
the radius of the WS wave function for both the s-wave 
and p-wave states in the limit of ews ~^ 0. In contrast, 
the HFB wave functions show only a small increase of 
radius even in the limit of ews ~^ 0. This feature remains 
the same even when the contribution of the other orbits 
are taken into account, as shown in the middle panel of 
Figs. 1 and 2. Due to the pairing effect in the continuum. 
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FIG. 1: (Color online) The mean square radii and the average 
paring gap as a function of the single particle energy ews in 
a Woods-Saxon mean-field potential. The top panel shows 
the mean square radius of the 2si/2 wave function with and 
without the pairing correlation, denoted by HFB and Woods- 
Saxon, respectively. The middle panel shows the rms radii 
for ^^O (the dotted line), ^^O (the dashed line), and ^''O (the 
solid line), obtained with the Hartree-Fock (^^O and ^^O) and 
the Hartree-Fock-Bogoliubov (^''O) calculations. The bottom 
panel shows the results of the HFB calculations for the average 
paring gap of ■^*0. 



the HFB calculations yield a larger radius than the HF 
calculations for the cases of ews < —1 MeV. On the 
other hand, in the case of —1 MeV < ews < MeV the 
HF wave function (equivalently one quasi-particle wave 
function in HFB) extends largely, while the HFB wave 
function docs not get much extension due to the pairing 
anti-halo effect. 

In the bottom panel of Figs. 1 and 2, the average pair- 
ing gaps are shown as a function of the single particle en- 
ergy ews • It is seen that the average pairing gap increases 
as ews approaches zero. See also Fig. 2 in Ref. [11] • This 
is due to the fact that the paring field couples with the 
extended wave functions of weakly bound nucleons in the 
self-consistent calculations. That is, the pairing field is 
extended as the wave functions do and becomes larger for 
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FIG. 2: (Color online) Same as Fig.l, but for the 2p3/2 state 
and for 30,3i,32]v^g isotopes. 



a loosely bound system. On the other hand, the paring 
gap behaves differently in the non-self-consistent model 
for the pairing field, in which the pairing field is fixed 
while solving the HFB equations [21|. In this model, 
the pairing field has a smaller overlap with the loosely 
bound nucleons because of the lack of the coupling, lead- 
ing to a smaller paring gap in the limit e — t- 0. Thus the 
anti-halo effect is somewhat underestimated in the non- 
self-consistent model as is expected from the asymptotic 
behaviour of the HFB wave function given by Eq. ([6]). 

Let us now calculate the reaction cross sections for the 
30,3i,32]\jg isotopes and discuss the role of pairing anti- 
halo effect. To this end, we use the Glauber theory, in 
which we adopt the prescription in Refs. [isl [2^ in or- 
der to take into account the effect beyond the optical 
limit approximation. Fig. 3 shows the reaction cross 
sections of the 30.3i.32]s^g nuclei on a ^^C target at 240 
MeV/nucleon. We use the target density given in Ref. 
[lOj and the profile function for the nucleon-nucleon scat- 
tering given in Ref. [l^. In order to evaluate the phase 
shift function, we use the two-dimensional Fourier trans- 
form technique [2^ . The cross sections aji shown in Fig. 
3 are calculated by using projectile densities constrained 
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FIG. 3: (Color online) Reaction cross sections of Ne isotopes 
on a ^■^C target at i5iab=240 MeV/A. The cross sections are 
calculated with the Glauber theory with HF and HFB densi- 
ties. The solid line with the filled squares shows the results of 
S'„(^^Ne) =|e(2p3/2)|=0.32 MeV, while the dashed line with 
the open triangles is obtained for 5'„('^^Ne) =|e(2p3/2)| = 1.49 
MeV. The experimental data are taken from Ref. [§]. 



to two different separation energies of the 2^3/2 neutron 
state. The dashed line with triangles is obtained using 
the wave functions with the separation energy Sn ('^^Ne) 
= jews I =1-49 MeV, while the solid line with squares is cal- 
culated with the wave functions of S'„('^^Ne)=|ews|=0.32 
MeV. The empirical separation energy of '^^Ne has a large 
ambiguity with S'„=0.29±1.64 MeV^. The cross sec- 
tion aR of ■^'^Ne is already much larger than the system- 
atic values of Ne isotopes with A<30. On top of that, we 
can see a clear odd-even staggering in the results with 
the smaller separation energy, S'„=0.32 MeV, as much 
as in the experimental data, while almost no stagger- 
ing is seen in the case of the larger separation energy, 
S'„=1.49 MeV. This difference is easily understood by 
looking at the anti-halo effect for |ews| ^ IMeV shown 
in the middle panel of Fig. 2. Recently, the effect of de- 
formation of neutron-rich Ne isotopes on reaction cross 
sections was evaluated using a deformed Woods-Saxon 
model [2^. It was shown that the deformation is large 
as much as (32 ^ 0.42 in '^^Ne and enhances the reaction 
cross section by about 5%. However, the calculated re- 
sults did not show any significant odd-even staggering in 
an between ^SNe and ^^Ne [26j . 

We mention that the pairing anti-halo effect may ap- 
pear also in lighter halo nuclei. In Borromean systems, 
only the three-body system as a whole is bound and none 
of the two-body subsystems is unbound. In the case of 
^^Li nucleus, this idea is implemented in that ^^Li is 
bound while ^''Li and a di-neutron are unbound. That is, 
the pairing correlations work in the continuum and gain 
the binding energy to make the three-body system bound 
[13, H^. We can make a simple estimate on how large the 
anti-halo effect is in the ^^Li nucleus. The experimental 
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FIG. 4: (Color online) Reaction cross sections of C isotopes 
on a ^^C target at _Eiab=83 MeV/A. The cross sections are 
calculated with the Glauber theory with three-body model 
densities. The experimental data are taken from Ref. [33 |. 



matter radii were obtained for ^Li and ^^Li as 



= 2.43 ± 0.02fm ^Li, 
3.27±0.24fm "Li, 



(8) 



respectively[l|. The matter radius of the two neutrons 
outside the ^Li core can be estimated by using a formula 
(see e.g., Ref. (2^]), 

(r^)2n = y ((r^)™("Li) - ^(r2)„(^Li)) = (5.68)^fm^ 

(9) 

On the other hand, the rms radii of single particle 
states lpi/2 and 2si/2 could be calculated with a Woods- 
Saxon potential adjusting the separation energy as Sn — 
|£wsh0.15 MeV; 



1^ = 6.48 fm for Ip^ 



/2 



'i^ = 10.89 fm for 2si/2 



(10) 



These values lead to the rms radius of ^^Li to be 3.53 
fm and 5.14 fm assuming pure (lpi/2)^ and (2si/2)^ two- 
neutron configurations outside the ^Li core, respectively. 
These values are apparently larger than the empirical 
value, 3.27 fm, and thus one sees the pairing anti-halo 
effect in the matter radius of ^^Li. 



The manifestation of the pairing anti-halo effect may 
be trivial for ^^Li as ^^hi is unbound. We therefore in- 
vestigate next neutron-rich C isotopes. We particularly 
study the ""^^C nucleus using a three-body model given in 
Ref. [lO]- In this case, the valence neutron in ^'^C oc- 
cupies the 2si/2 level at ews = —1.21 MeV, while ^^C 
is an admixture of mainly the (2si/2)^ and (1^5/2)^ con- 
figurations. Assuming the set D given in Ref. )30,] for 
the parameters of the Woods-Saxon and the density dis- 
tribution for ^^C given in Ref. [3l|, the rms radii are 
estimated to be 2.53, 2.90, and 2.81 fm for "C, ^^C, and 
^^C, respectively. The corresponding reaction cross sec- 
tions an calculated with the Glauber theory are shown in 
Fig. 4. The calculation well reproduces the experimental 
odd-even staggering of the reaction cross sections, that 
is a clear manifestation of the pairing anti-halo effect. 

In summary, we have studied the mass radii of Ne 
isotopes with the Hartree-Fock (HF) and Hartree-Fock- 
Bogoliubov (HFB) methods with a Woods-Saxon poten- 
tial. The reaction cross sections were calculated using 
the Glauber theory with these microscopic densities. We 
have shown that the empirical odd-even staggering in the 
reaction cross sections of neutron-rich Ne isotopes with 
the mass A=30 ~ 32 is well described by the HFB den- 
sity and can be considered as a clear manifestation of the 
pairing anti-halo effect associated with a loosely-bound 
2p3/2 wave function. The anti-halo effect was examined 
by the three-body model calculations also for ^^Li and 
^^C. We argued that the observed mass radius of ^^Li 
can be considered as an evidence for the anti-halo effect. 
We have also shown that the experimental reaction cross 
sections of i4-16q a,re well reproduced, confirming the 
ani-halo effect in these nuclei. 

We point out that the odd-even staggering in reaction 
cross sections has been widely observed in a series of nu- 
clei with small separation energies, such as 14-16(^ [33 |. 
i8-20ci, 28-30Ne, 30-32Neg and 36-38Mg[3|. Evi- 
dently, the pairing anti-halo effect is commonly seen in 
those weakly-bound nuclei near the drip line. 
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